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Abstract. We report first results from an ongoing program of combining visibilities from the Giant Meterwave Radio Telescope 
(GMRT) and the Nancay Radio Heliograph {NRH) to produce composite snapshot images of the sun at meter wavelengths. We 
describe the data processing, including a specific multi-scale CLEAN algorithm. We present results of a) simulations for two 
models of the sun at 327 MHz, with differing complexity b) observations of a complex noise storm on the sun at 327 MHz on 
Aug 27 2002. Our results illustrate the capacity of this method to produce high dynamic range snapshot images when the solar 
corona has structures with scales ranging from the image resolution of 49 to the size of the whole sun. 
We find that we cannot obtain reliable snapshot images for complex objects when the visibilities are sparsely sampled. 
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1. Introduction 

We have an ongoing program of combining data from the 
GMRT and NRH in order to make composite meter wavelength 
images of the sun. Our aim is to combine the short baseline 
data from the NRH with long baseline data from the GMRT in 
order to obtain meter wavelength images of unprecedented res- 
olution and fidelity. We present here the first results from our 
campaign. 

Imaging complex, rapidly varying radio sources in the solar 
corona is problematic because one needs to obtain images span- 
ning durations as short as a few seconds. The uv-coverage for 
short duration exposures is usually quite poor, and this affects 
the ability to image complex sources reliably. Furthermore, de- 
convolution procedures like the commonly used CLEAN often 
perform poorly on sources that exhibit structures on a variety 
of scales (Wakker & Schwarz 1991). The problem of poor uv- 
coverage can be alleviated by synthesis imaging over time du- 
rations of a few hours, but this technique cannot be used for 
studying phenomena that vary rapidly in time. There are sev- 
eral instances in solar physics where high dynamic range, high 
resolution snapshot images are essential. For instance, stud- 
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ies of extended thermal emission in the corona, coronal mass 
ejections (CMEs) that frequently occur together with compact, 
bright bursts and flare continua require imaging with a large 
field of view, high spatial resolution and dynamic range. High 
resolution observations can establish a firm limit on the angu- 
lar broadening of sources in the solar corona due to turbulence, 
which determines the apparent size of radio bursts (e.g. Bastian 
1994). 

1.1. Dynamic Range 

Before proceeding further, we pause to define two kinds of dy- 
namic ranges we will refer to below. The rms dynamic range 
is the definition usually employed by the community. It is the 
ratio of the peak brightness of the image to the rms brightness 
of a representative empty region near the area where the peak 
brightness occurs. This definition is relevant when artefacts on 
the final images look like random noise, which is often the case 
for images produced via aperture synthesis. However, we deal 
here with snapshot images, for which the visibilities are rather 
sparsely sampled, especially for large baselines. The artefacts 
on such snapshot images are localized and structured, and of- 
ten look like spurious compact sources. We will use another 
definition of dynamic range, which we call the "max" dynamic 
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range, which we think is more relevant to the snapshot images 
we present. The max dynamic range is defined as num/denom, 
where num. is the peak brightness of the image. The quantity 
denom is the maximum deviation from the mean of an empty 
region around the area where the peak brightness occurs. We 
will quote the rms dynamic range for the sake of comparison 
with pubUshed results. The max dynamic range is typically an 
order of magnitude lower than the rms dynamic range for our 
results. 

1.2. Meter wavelength observations of the solar 
corona: a brief overview 

Meter wavelength observations of the solar corona to date have 
concentrated either on having a large field of view with dense 
short baseline coverage (e.g., observations using the NRH) or 
on achieving high spatial resolution at the expense of a sparse 
short baseline coverage (e.g., extended array observations with 
the Very Large Array (VLA)). In what follows, we attempt to 
give a short overview of meter wavelength imaging observa- 
tions of the solar corona. This will help to place the observa- 
tions reported here in context. 

There have been several attempts at detecting small scale 
structures at meter wavelengths in the solar corona. Kerdraon 
(1979) used one-dimensional (ID) scans at a time cadence of 
0.02s and a spatial resolution of 1 .2' with the NRH at 1 69 MHz. 
The ID scans were obtained by a simple Fourier transform 
(FT) of the observed complex visibihties (CVs). The field of 
view was around 1 degree, and the artefacts on the ID scans 
were due to sidelobes of the response function. The "max" 
dynamic ranges of these ID scans were 5-10. Kerdraon con- 
cluded that type 1 bursts at 169 MHz range in size from 0.7' to 
3'. 

Zlobec et al (1992) have used the VLA in its A configura- 
tion with a resolution of 4 " at 327 MHz to search for small scale 
structures in snapshot maps with a time cadence of 1.6 seconds. 
They find that the smallest type I source sizes they can observe 
are around 30 ". The largest structure that can possibly be im- 
aged with the VLA A- array is around 110 . However, owing 
to the shortage of short baseUne spacings, Zlobec et al. (1992) 
mention that they cannot be confident about large scale struc- 
tures. In some cases, they cannot be confident about the po- 
sitions of the small structures. Furthermore, the rms dynamic 
range of the images they obtain is around 10, at best. 

Some other instances of snapshot imaging of noise storm 
sources in the solar corona at meter wavelengths are Willson 
et al. (1997); (VLA C configuration, resolution of 40 " at 327 
MHz, time cadence 3.3 seconds), Willson et al. (1998); (VLA 
A configuration, resolution 30 " at 75 MHz, 5 " at 327 MHz, 
time cadence 10 seconds, source size 1.5' at 327 MHz, 2.5 x 
5' at 75 MHz), Willson (2000) (VLA D configuration, time ca- 
dence 3.3 seconds, resolution 180 at 327 MHz). 

There have also been several instances of synthesis imaging 
of noise storm emitting regions at meter wavelengths. Habbal et 
al. (1996) have used the VLA in its B configuration to observe 
noise storm sources with synthesis durations ranging from 
around 100 mins to 3.5 hours. Their resolution was around 94 " 



at 327 MHz. EarUer, Habbal et al. (1989) had used the VLA in 
its C configuration to image a noise storm source at 327 MHz 
with a synthesis duration of around 4 hours. They were able to 
image structures ranging from around 80 ' to 240 ' with a reso- 
lution of 57" X 47" at 327 MHz. Willson et al. (1997) have used 
data from the VLA C array to obtain 4 hour synthesis images 
of a faint transequatorial arch measuring 600" at 327 MHz. In 
a multiwavelength study aimed at high dynamic range imag- 
ing. White et al. (1992) have used the VLA D configuration to 
obtain 4 hour synthesis images of active region plages measur- 
ing around 600 with a resolution of around 200 at 327 MHz. 
The rms dynamic range of their 327 MHz images was between 
60 and 90. In a similar study at 1420 MHz, Gopalswamy et 
al. (1991) used the VLA D configuration to obtain 6-7 hour 
synthesis images of active region plages ranging from 50 to 
200" with a resolution of 50". The rms dynamic range of their 
images varies from 170 to 346. 

1.3. Present Study 

In this paper, we present composite snapshot images by com- 
bining visibilities from the GMRT and the NRH for : 

- simulations with two models of the radio sun at 327 MHz. 
The first model is relatively simple and similar in many as- 
pects to the real case reported in this paper. The second 
model is more complex, and we use this in order to investi- 
gate the limits of what we can hope for with snapshot imag- 
ing. 

- observations of complex and evolving noise storm emitting 
regions in the solar corona at 327 MHz on Aug 27 2002. 
These observations use an integration time of 17 seconds. 
They illustrate the capabilities of the composite instrument 
with real data, but since solar bursts vary over time scales 
that are much shorter than 17 sec, they are of limited use 
in arriving at astrophysical results. 

This procedure complements the dense short basehne coverage 
of the NRH with the long baseUne coverage of the GMRT. We 
obtain snapshot maps of structures with sizes ranging between 
the resolution (49 ) and the size of the whole sun with unprece- 
dented (rms) dynamic ranges of 250^20. 

2. The Instruments 

The Nancay RadioHeliograph (NRH) and the Giant Meterwave 
Radio Telescope (GMRT) both operate at meter wavelengths, 
and we currently have common observations of the solar corona 
at 236 and 327 MHz. The two instruments have up to 4 hours of 
common observing time each day. The following subsections 
give brief overviews of the two instruments. 

2.1. The Nancay Radio Heiiograph (NRH) 

The NRH (Kerdraon & Delouis 1997) is located about 200 km 

south of Paris (latitude 47''23'A^, longitude 2° 12' E, altitude 
150 m). It consists of 44 antennas of size ranging from 2-10 
m, spread over two arms {EW and NS) with respective lengths 
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of 3200 m and 2440 m. The EW and NS baselines dsw and 
Jjvs are from to 3200 m (resp 2440 m) by steps 100 m (resp 
54 m). After January 2003 the number of antennas were raised 
to 48 and the EW minimum baseline is reduced to 50 m. The 
NRH observes the sun 7 h per day at up to 10 frequencies be- 
tween 150 and 450 MHz. The standard observing mode con- 
sists of Stokes I and V images at 164, 236, 327, 410 and 432 
MHz each 150 msec. The resolution of two dimensional (2D) 
images and the field of view depend on the frequency and the 
season. For the observation reported here they are ~ 4.8 x 3.0 
and ~ 40 x 70 respectively. 

2.2. The Giant Metrewave Radio Telescope (GMRT) 

The GMRT is located about 80 km north of the city of Pune 
in Maharashtra, western India (latitude 19"06 N, longitude 
74"03 E, altitude 650 m). It consists of thirty 45 m diameter 
antennas spread over 25 km. Half of these are in a compact, 
randomly distributed array of about 1 square km and the rest 
are spread out in an approximate ' Y' configuration. The GMRT 
currently operates in the frequency bands around 235, 325, 610 
and 1000-1450 MHz, with bandwidths ranging from 16 to 32 
MHz. Of these, the 235 and 327 MHz frequencies are most 
commonly used for solar observations. The shortest baseline is 
100 m and the longest one is 26 km. Data is usually acquired 
at a time cadence of 16 seconds, but it is presently possible to 
acquire data at intervals as short as 2 seconds. The automatic 
level control circuitry in the GMRT data processing electronics 
has a response time of around 1 second. This limits the short- 
est integration time to be 2sec, which is still substantially 
longer than what can be achieved with the NRH. Solar ob- 
servations are typically caiTied out from around 03:30 UT to 
12:30 UT. During such observations, the sun is typically ob- 
served for around 15-20 minutes, after which a nearby phase 
calibrator is observed for around 10 minutes. This cycle is then 
repeated. During a particular cycle, the antennas track the sun 
at a linear rate in RA and DEC. At the beginning of the next 
cycle, the antennas are repositioned to the RA and DEC of the 
sun corresponding to the time when the cycle begins, and they 
then track the sun as before. Further details about the GMRT 
can be found at the URL http://www.gmrt.ncra. tifr.res.in|and in 
Ananthakrishnan & Rao (2002), Swarup et al. (1991).^ 

2.3. Tiie composite instrument 

It is useful to discuss at this point some properties of the com- 
posite instrument, for a better understanding of the subsequent 
data processing. Fig. 1 shows the uv coverage of the joint in- 
strument for 2002 aug 27 at 09:30 UT, with baselines up to 12 
000.1: 



+++ 




(a) 



i_LL 



1. The NRH produces the central regular pattern (see fig. lb): 
a densely covered parallellogram with some linear exten- 

smaller than 



sions, limited to values of r,, 



=s 3000/1. This results from the regular antenna distribution 
oflhsNRH. 

The GMRT produces an irregular and much more extended 
pattern, consisting of: 



(b) 



Fig. 1. Actual iiv coverage (rac/-') for 2002 Aug 27, 09:30 UT 
with NRH + GMRT. a) coverage with basehnes up to 12 000 
A, b) zoom-in of the central area showing the overlap between 
NRH and GMRT coverages. GMRT data are marked by +, 
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- a central set of points corresponding to baselines from 
the central compact array of GMRT, which has a sub- 
stantial overlap with the parallellogram of the NRH. We 
refer to the area of overlap between the inner baselines 
of the GMRT and NRH as ^overlap- 

- similar patterns with fewer points, distributed over 6 ir- 
regular arms,which correspond to baselines between the 
central array and the arm antennas. 

- isolated points corresponding to baseUnes between the 
arm antennas. 

The composite dirty lobe is the FT of this uv coverage and is 
the sum of NRH and GMRT dirty lobes rh and L^gmrt ■ Due 
to the regular NRH uv coverage Lcirh comprises of a pencil 
beam with classical "sine" side lobes repeated over a period 
« 1°. (larger than the size of the sun) on both axes (Fig. 2a). 

Within one period, most of the flux of Lj rh is concentrated 
in a small area surrounding the main pencil beam, with extent 
« 2 or 3 times the width of the pencil. The resulting dirty image 
of the sun looks like a smoothed version of the sun, with side 
lobes < 20% near compact sources. 

The situation is somewhat different for LdOMRT- The 
GMRT uv coverage is much more sparse and irregular and the 
beam is very complex (see fig. 2b) : the narrow main pencil 
beam contains only a small fraction of the flux and is sour- 
rounded by a complicated pattern of large side lobes. At larger 
distances from the main pencil beam, the structure of lobe is 
aperiodic (no regular grating lobes as for the NRH) and ex- 
hibits features of various scales and amplitudes, the most in- 
tense ones being at distances somewhat smaller than the size of 
the sun. As a result, the GMRT alone : 




(a) 



- cannot give a reliable snapshot image of the whole sun, es- 
pecially of the large scale features, 

- the dirty image of several sources on the sun is complex 
and difficult to understand visually, and requires a cleaning 
procedure. 

Similarly, the benefits of merging GMRT data with NRH data 
is not clear from a visual inspection of the dirty composite im- 
ages (which appear visually to be similar to NRH dirty images 
plus some "small scale noise"). The benefits are obvious only 
from the composite clean images. 



3. Data selection and processing 

We present observations of the sun at 327 MHz on Aug 27 
2002 between 9 and 10 UT in this paper. We chose this day 
because there are complex sources on the sun, which present a 
challenge for snapshot imaging. The GMRT observations were 
integrated for a duration of 17 seconds, while the NRH data 
were integrated for a duration of 0.125 seconds. 



3.1. Pre-processing 

The GMRT data were separately pre-processed in the follow- 
ing manner, using the standard AIPS software: a suitably inter- 
ference free channel was first identified from the 128 channels 



(b) 



Fig. 2. Theoretical lobes for (a) NRH and (b) GMRT. The cir- 
cle represents the solar optical limb. 
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spanning 16 MHz around 327 MHz. The data were then rig- 
orously edited for bad antennas and baselines. The visibilities 
on the sun were then caUbrated using those from the nearby 
phase calibrator 1021+219. Before any subsequent processing, 
the NRH data were therefore integrated over the 17 second in- 
tervals defined by the GMRT data. The high time resolution 
NRH data were first integrated over the 17 second intervals de- 
fined by the GMRT data, in order to match the time cadence 
of the GMRT data. The NRH and GMRT visibihties were then 
intercalibrated in the following manner: 

3.1 .1 . Amplitude intercalibration 

The visibilty amplitudes from NRH and GMRT should be co- 
herent, since both instruments use cosmic calibrators. NRH 
and GMRT uv coverages overlap particularly densely on the 
area Aoveriap around the origin of the uv plane (corresponding 
to baseUnes < Ikm). Hence the sets of samples of the visibil- 
ity on this area from both instruments should have comparable 
amplitude and merge homogeneously. This is especially true 
for the sun, for it dominates the system temperature for both 
the individual instruments. We noted that the amplitudes of the 
GMRT visibilities were larger (as compared to the NRH ones) 
by roughly one order of magnitude. This is primarily because 
the GMRT antennas are not fitted with solar noise calibrators, 
and the so-called "T^y" correction (Bastian 1989) can therefore 
not be appUed. In other words, we cannot be certain about the 
absolute values of amplitudes of the GMRT visibilities (when 
observing the sun). We applied a multiplicative factor aintei-cai 
to the GMRT visibilities. The value of aintereai was determined 
in such a way that the mean value for the modulus of visibility, 
when averaged over Aoveriap, was the same for both sets of data. 
We estimate that the accuracy of the resulting intercalibration 
is a few %. Visual inspection of amplitude versus uv distance 
plots for the composite data convinced us that this simple pro- 
cedure gives satisfactory results. 

3.1.2. Position intercalibration 

The absolute positions of sources as observed by the two in- 
struments can difl'er for several reasons: 

1. For merging the visibilities, data are extracted from the 
standard processing algorithms of each instrument. Since 
fringe stopping for solar observations is complex and not 
necessarily done in one step, this can produce slowly vary- 
ing position diff'erences between NRH and GMRT data at 
this stage. 

2. The inhomogeneities due to gravity waves in ionospheric 
electron density can produce apparent shifts of the order 
of l' at 327 MHz, varying on time scales ~ 30 min (see 
for instance Mercier, 1996 and references therein). These 
eff'ects are of local origin and would be totally uncorrelated 
between NRH and GMRT sites. 

It is obvious that the position intercalibration must be better 
than the expected resolution. This was achieved in the fol- 
lowing manner: for each 17 second interval, 2D dirty images 
were separately calcultated for NRH and GMRT through a 
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fast Fourier Transform (fl^t), including only baselines in the 
range common to both instruments (baselines < 2 km, which 
makes this range somewhat larger than Aoveriap). in order to 
achieve comparable resolutions. We then defined the transla- 
tion required to make the maximum of GMRT dirty image co- 
incide with the maximum of NRH dirty image. This enables us 
to calulate the phase shifts that must be applied to all GMRT 
baselines in order to make the spatial origin identical for both 
instruments. We find that these differences in position were ^ 
0.15 Rs and vary by less than 0.02 Rs for the duration of our 
observations (9-10 UT, Aug 27 2002). 

3.2. Processing 

The processing comprises essentially of two steps : getting a 
"dirty" image and cleaning it, in order to make it visually un- 
derstandable. When a single intense compact source is present, 
we can also estimate its dimensions by fitting a 2-D Gaussian 
directly to the observed visibilities. As we will see below, this 
allows us to measure sizes smaller than what would have been 
possible via inspection of clean composite images. 

3.2.1. Dirty image 

The dirty image is merely the Fourier Transform (FT) of the 
joint set of visibilities. The ff't algorithm has the advantage of 
being efficient, but it requires that spatial frequencies corre- 
sponding to the various baselines are distributed on a regular 
square grid in the uv plane. This is not the case because NRH 
and GMRT basehnes are not simply commensurable (GMRT 
antennas are randomly distributed). This situation is usual in 
radioastronomy and regridding is commonly used. It consists 
of replacing the observed samples of the visibility by samples 
interpolated on a regular grid. We show hereafter that it is not 
convenient for snapshot images of a complex and wide object 
such as the sun with NRH and GMRT. 

According to Shannon's theorem, if w is the total angu- 
lar width of the sun, its visibilty must be sampled with a step 
Usampie which Satisfies : 

^sample ^ (1) 
W 

Taking w ~ 45' = 0.013 rod at 327 MHz, Usampie = llrad'^. 
The corresponding distance between neighboring antennas 
(projected normal to the line of sight) is i x UsampU - 69m. 
The sampling condition (1) is fuUfiUed for NRH baselines (at 
least far from local noon for observations before January 2003) 
but not for the GMRT baselines, especially for large baselines 
involving the distant arm antennas (which yield the enhanced 
spatial resolution). Thus interpolated values of the visibility for 
large baselines can be strongly biased and produce spurious ef- 
fects such as "ghosts" of the actual sources. We emphasize that 
the condition (1) should be strictly satisfied for interpolating 
the visibility, since there can be several sources on the sun that 
are separated by more than 2 Rs . 

In order to overcome this difficulty, we have defined a grid 
in the uv plane with a period Ugrid « Usampie- In this way the 
observed visibilities can be ascribed to the nearest node of the 
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grid, with only a small error. In other words, we keep the ex- 
act values measured for the visibility and make a small error 
("round-up" error) < Q.Sugrid on the spatial frequencies, which 
is much smaller than the minimum possible scale of variation 
of the visibility in the uv plane, according to (1). In this way, 
we avoid the problem of interpolation in the uv plane, at the 
expense of handling large arrays in fft computations. 

In our data processing we used arrays of dimension 1536 x 
1536 and restricted the longest baseline to 17000/1 (beyond 
this value GMRT uv coverage is too sparse for snapshots, see 
§ 3.2.2 below). Hence Mg„v/ ~ llrad^^ and the round-up er- 
ror in spatial frequencies is < llrad^\ which is « Usampie- 
Since NRH and GMRT taken together usually give fewer than 
1000 visibility samples, most of the elements of the 1536 x 
1536 array in the uv plane are zero. When computing the im- 
age through a fft, the resulting field is l/ugrid ~ ^deg.. This 
field is much larger than the maximum width of the sun observ- 
able with the NRH (see eq.(I)), and we therefore retain only its 
central part. 

3.2.2. Cleaning procedure 

We first briefly recall the basic idea of CLEAN deconvolution. 
When the uv coverage is sparse and irregular, the dirty lobe Ld 
has a compUcated shape, with large side lobes at distances from 
the main peak that are much larger than its half power beam 
width (HPBW). The contribution of a point on the object is not 
locaUzed within one HPBW of Ld on the dirty image, and this 
makes the dirty image difficult to interpret visually. CLEAN 
replaces the convolution of each point with the dirty lobe by a 
convolution with a clean lobe (which is usually a Gaussian, 
with no sidelobes). The amplitude and HPBW of the clean lobe 
are similar to that of the dirty lobe. The procedure is iterative, 
starting with the brightest points on the dirty image. The clean 
image thus depends "locally" (i.e., within one HPBW) on the 
object and is thus visually more understandable. For technical 
reasons the replacement is only partial at each step (see for in- 
stance Thompson et al, 1986, chapter 11). In the uv domain, 
the CLEAN process can be viewed as one where the holes in 
the uv coverage are filled up via interpolation. However it is 
well known that CLEAN works rather poorly for continuous 
objects that have a wide spectrum of spatial scales. Some reme- 
dies have been proposed to alleviate this problem. In particular, 
when the object involves two well separated spatial scales (eg., 
few compact sources superposed on the quiet sun), a two-scale 
CLEAN procedure can be used (Wakker and Schwarz 1991). 
One can also try to subtract a model of the quiet sun from the 
dirty image. CLEAN then works on a dirty image which sup- 
posedly results only from only some isolated compact sources. 
The model of quiet sun is added afterwards to the cleaned im- 
age for obtaining a final clean image. 

The cleaning procedure we developed here is a multi-scale 
version of CLEAN (MS - CLEAN) adapted to objects with 
features of various scales (Delouis 1998). The basic idea is to 
first produce several filtered dirty images Ido, ... Idn from the 
original dirty image /. Each filtered dirty image involves a rel- 
atively narrow range of spatial scales. For this purpose, we use 



successive filters Fq, F„ in the uv plane, each filter selecting 
progressively smaller spatial scales: 

- Fo is a low pass filter with cut-ofl" Ueo ~ 3^ Usampie (recall 
that Usampie ~ 3.5 Ugrid and that the observable field of view 

is 1 / U sample) • 

- The filters F\, ...F„ filters are band pass filters with internal 
radii un - Ueo, ... Mm = fl"~'MeO and external radii u^i - 
aua = aueo, ■■■ Ue„ = a"ui, where a 3 in order to limit the 
spatial scale ranges. 

The edges of these filters are smooth and their sum is unity. 

The first filtered dirty image Ido is very smooth and com- 
prises of only a few independent pixels. It can be cleaned using 
smooth dirty and clean beams Lda and L^o deduced from Ld and 
Lc through a filtering with Fq. Thus the largest structures in the 
dirty image can be described in terms of only a small number 
of independent smooth components. The problem of describing 
a large scale structure with a large number of narrow compo- 
ments (if the CLEAN algorithm were directly applied) is then 
reduced to that of describing a source in components « 4 times 
narrower than the widest source that can be described. It is ob- 
vious that the number of iterations (and with it, the possibility 
of producing artefacts) in the CLEAN procedure will be dras- 
tically reduced. 

The same argument can be applied to other filtered dirty 
images : each case is not very different from the case of a rela- 
tively small number of compact sources ("compact" relative to 
the dirty lobe). The final clean image is simply the sum of the 
filtered clean images. 

After some trials we found that the best results were ob- 
tained under the following conditions: 

- The GMRT uv coverage gets increasingly sparse for pro- 
gressively longer baselines. At 327 MHz there are only 
isolated points for baselines > 17000/1. This could lead 
to artefacts when cleaning dirty images from complex ob- 
jects. We found that the best compromise between reducing 
these artefacts and preserving the spatial resolution could 
be achieved by taking a clean lobe Lc for which the width 
of its FT is » 0.35 X uVmax = 0.35 X bmaxl^, where b^ax is 
the longest baseline under consideration. In simulations we 
could take uVmax - 17000 or 12000, whereas for real data 
we had to take uv^ix = 12000 only, probably because of 
imperfect calibration. In other words, MS - CLEAN works 
poorly for sparse uv coverages. The resulting HPBW of the 
clean lobe was » 1 . 

- We took a maximum index n=3 for the filters (i.e., a central 
zone and 3 concentric rings). This allows us to cover a spa- 
tial frequency range up to » IQ'^rad'^ with a ratio a » 3 for 
the ratio between the external and internal radii of filters. 

3.2.3. Estimating the width of an intense compact 
source 

The resolution corresponding to the width of our final clean 
beam is much lower than what could be expected from the 
GMRT size. As explained earlier, this limitation arises because 
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we try to obtain images from a complex object with a cover- 
age that becomes very sparse for long baselines. However, for 
particularly simple objects, (e.g., when the composite image 
contains a single intense, compact source) it is possible to esti- 
mate the source size by fitting a gaussian model directly to the 
amphtude of the visibihty over the whole uv range observed 
by the GMRT. It is then possible to compare this derived size 
with what is directly measured on clean images. Of course, this 
method is meaningful only when there is a dominant compact 
source, and cannot be appUed when several sources with simi- 
lar intensities co-exist. 

4. Results 

We first checked the expected performance of our method us- 
ing simulated brightness distributions that are similar to the ac- 
tual distribution on the day of our joint observations. We found 
that this procedure gives us a reliable method of anticipating 
the limitations when we process actual data. We show here the 
results of two representative simulations. 

4.1. Simulations 

4.1 .1 . The models for the radio sun 

We defined a class of solar models which contain several fea- 
tures commonly observed at 327 MHz. They comprise struc- 
tures of various sizes, ranging from the size of the quiet sun, 
down to 0.05 Rs - 52 . These structures are derived from a 
simple mathematical model. The brightness distribution in the 
ID analog of this model can be written as 

bix) ~ 1^ ~ |~ j j ' -a < X < a , 

bix) = 0, \x\>a. (2) 

The parameters a, ei and 62 define the half-power width of 
the structure, the flatness of its maximum and the slopes of its 
sides, respectively. In practice, we took a 2D model with an 
elliptical section, where the tilt of the ellipse major axis is also 
a free parameter. The model of the sun is the sum of ^ 30 
such structures with various positions, sizes, flatness, tilt and 
amplitudes. The structures are: 

- a relatively flat structure for the quiet sun, with total extent 
42 x38' (EW X NS), and Tb = lO^K. 

- several "thermal" regions with various slopes, widths in the 
range 6 - 9 , various tilts and excess brightness of some 
10^ A" above the quiet sun. 

- two coronal holes, one very elongated, with steep sides and 
depression depths 3 x 10^ K, 

- a number of compact elliptical "non thermal" sources with 
sizes 1-3', various tilts and excess brightness of some 
lO^K-3 X 10'^ K above the quiet sun, 

- a faint CME-like arch whose brightness decreases with 
height and with distance from the central axis. 

Using these ingredients, we produced two particular mod- 
els for our simulations : 




Fig. 3. Model 1 adopted for the simulations. 

1. A relatively simple model (model 1) in which only 4 "non 
thermal" sources are intense (T up to 3 x 10^ K) and there- 
fore dominate the other structures, 

2. A more complex model (model 2), which has "non ther- 
mal" sources with T ~ lO^K - 5 X lO^K, so fliat they do 
not dominate the other features as in model 1. 

4.1.2. Model 1 

Fig. 3 shows model 1. It is similar to several instances encoun- 
tered in the real data of Aug 27 2002. It consists of : 

- a bright western source with Tb = 3x lO^K, 

- two relatively dimmer sources near the disk center and 
extended emission (of size 4 ) between these two central 
sources. 

[The color scheme for fig. 3 (and the subsequent figures) was 
chosen so as to optimize its visual appearance on paper prints. 
The intensity scale is indicated on the lower left corner: it is lin- 
ear for low levels and progressively saturates at higher levels. 
This highlights low level features {CLEAN artefacts, in partic- 
ular)]. 

All other features (including the quiet sun) are hardly visi- 
ble. Fig. 4 shows the clean images obtained by the NRH alone 
(a), GMRT alone (b) and NRH + GMRT (c) when using the 
actual uv coverages provided by the NRH and GMRT at 09:30 
UT on Aug 27 2002. This means that the GMRT baselines 
which were flagged for the actual data of Aug 27 2002 were 
not used here either. The synthesized clean beam is shown at 
the bottom left comer of each figure. For cases b and c, it has 
an extent of 49 . Even for this relatively simple model, the im- 
provement in quality of the composite clean image (c) is quite 
evident : the resolution is much better than the NRH-only im- 
age, and there are fewer spurious large scale structures than in 
the GMRT-only image. 

4.1.3. Model 2 

Model 2 is shown in Fig. 5. It differs from model 1 in that the 
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(b) 



(b) 




(c) 



Fig. 5. a) and b): Model 2 adopted for the simulations: a more 

complex sun. Panel a) shows an overview of the model. The 
color scheme has been chosen to enhance lower levels in panel 
b). 



the nonthermal sources are dimmer, which allows us to discern 
the other features better: the quiet sun, patches of "thermal" 
emission, coronal holes, faint compact sources and a CME-like 
arch in the SW. 

Fig. 6 shows the clean images obtained by the NRH alone 
(a), GMRT alone (b) and NRH + GMRT (c), using the same 
uv coverage as for model 1 : 

The NRH image is smooth (faint compact sources are not 
seen), but large scale structures are reliable, although they 
appear with a reduced contrast. The faint CME beyond the 
south west limb is satisfactorily evident. Conversely, due to the 
greater complexity of the object, the clean image from GMRT 
alone is practically unusable : large scale structures cannot be 



Mercier et al.: Combining visibilities from GMRT and NRH: High dynamic range 327 MHz snapshots 



9 





(a) 



Fig. 7. Clean composite {NRH + GMRT) image for model 2, 
using all the 435 possible GMRT baselines 




(b) 




(c) 



recognized, and only the two most intense compact sources can 
be identified. The improvement in quality of the composite im- 
age (c) is also evident, although hmitations due to the sparse uv 
coverage more apparent than for model 1 : 

- the edges of coronal holes are blurred 

- low level artefacts (both positive and negative) appear 
around coronal holes and thermal regions, 

- some of the faintest compact sources are difficult to distin- 
guish from artefacts. 

The CME is less visible on fig 6c than on fig 6a, but this 
is only due to the higher contrast of fig. 6c, and because of its 
superior resolution. 

Fig 7 shows the composite clean image obtained with the 
uv coverage that could be provided if all GMRT antennas were 
functioning. It is somewhat better than fig. 6c : faint compact 
sources are more reliable and the northern coronal hole has 
fewer artefacts. 

As expected, the benefits of merging NRH and GMRT 
data are obvious when compact sources are present. However, 
it may be noted that low contrast extended features (e.g., the 
coronal holes of model 2) are also better imaged with the com- 
posite instrument. Model 2 illustrates the hmits of what can 
be hoped from composite snapshots of a complex object with 
sparse extended uv coverage. 

In the preceding simulated examples, we had introduced 
no calibration errors (i.e., deviations in gains and phases). We 
have also produced images including random phase deviations 
in the simulated visibilities as a simple means of approximating 
the effect of phase calibration errors. We found that calibration 
errors are the principal cause of clean artefacts on the final im- 
ages, and can effectively limit the resolution : for nns phase 
deviations > 20 deg, clean artefacts become noticeable and the 
extent of the accepted uv coverage needs to be reduced. 
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4.2. Observed data 

We now turn our attention to the data actually observed be- 
tween 09:02:23 UT and 09:57:54 UT on Aug 27 2002. During 
this time, there was a strong ongoing noise storm near the west 
limb (0.69 Rs west, 0.34 Rs south), hereafter referred to as 
"W. There were also two weaker noise storm sources near the 
disk center (0.13 Rs east, 0.12 Rs south and 0.13 Rs east, 0.34 
Rs south), referred to as "A^' and "5 " respectively. While the 
western noise storm source remains compact throughout the 
duration of our combined observations, the two sources near 
disk center evolve significantly, showing marked extensions 
at certain times. We have selected three 17 second snapshots, 
which highhght some of these features and illustrate the capa- 
biUty of our combined imaging technique. 

Figure 8 shows 17 second snapshots centered around 
09:20:25 UT, using data from NRH alone (a), GMRT alone 
(b) and from NRH + GMRT (c). The uv coverage has already 
been shown in fig. 1 : the relatively sparse, extended uv cover- 
age of the GMRT, which is complemented by the dense, short 
baseline uv coverage of the NRH. 

The improvement in resolution (when the GMRT data are 
used) is evident : the half-power size of the clean beam in fig. 
8b and 8c (lower left hand comer) is 49 " . The rms dynamic 
range for the composite image (fig 8c) is 244, while the max 
dynamic range is 11. On the GMRT -only image (fig 8b), we 
note that : 

- artefacts are spread over the whole field. The relatively 
strong level of negative artefacts produces a grey back- 
ground for the picture. 

- the two central sources are not properly imaged. This is 
especially true of "N" which is smooth in fig 8b, whereas it 
is more compact on the composite image (fig 8c). 

The composite image shows artefacts which are both fainter 
(as is evident from the darker background) and more localized 
to the vicinity of sources. This must be compared to simula- 
tions with model 1 (which is very similar to this actual case), 
for which GMRT alone gives satisfactory results for compact 
sources and a much lower level for artefacts (fig. 4b). Since the 
data processing is the same for the simulations and for real data, 
we think that the difference should be ascribed to imperfect 
phase calibration of GMRT data, and possibly to phase shifts 
arising from the 30 dB solar attenuators, which are switched on 
for solar observations and off for calibrator observations. 

The western source "W" is brightest at 09:13:39. Fig. 9 
shows clean snapshot images at this time from NRH alone (a), 
GMRT alone (b) and from NRH + GMRT (c). Clearly, the im- 
proved resolution is provided by the GMRT baseUnes, while 
NRH data helps in reducing clean artefacts. Note that, as for the 
preceding case, 'W is also poorly imaged by the GMRT alone. 
The resolution of the composite image, as in earlier cases, is 
49 . The rms dynamic range of this image is 370 and the max 
dynamic range is 21 . 

The source "W" in this snapshot is the smallest structure 
present in our observations. It is barely resolved with the clean 
beam size of 49 . However, since there is only one intense 
source present at this time, we can use the step described in 
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§ 3.2.3 (fitting a gaussian to the visibilities) to estimate a size 
of a: 1 . Using simulations, we checked that this derived size 
is not appreciably affected by the faint sources 'W" and "5" 
and by the faint halo surrounding "N". High time resolution 
NRH data show that the sources 'W and "5" (and "IV", to a 
lesser extent) exhibit continual small-scale motions. Since the 
time resolution of composite images is limited by the 17 sec- 
ond integration time of the GMRT data, it is possible that this 
snapshot image has "smeared over" such small scale motions. 
It follows that the source size of I we derive for "W" is likely 
an over-estimate of its actual size. 

Fig. 10 shows 17 second snapshots centered around 
09:04:04 on Aug 27 2002, from NRH alone (a), GMRT alone 
(b) and from NRH + GMRT (c). It illustrates the abihty of 
the composite instrument to resolve and image extended emis- 
sion reliably. The extended emission of around 200" that lies in 
between, and to the west 'W and "S " is imaged much better 
in the composite image than in any of the individual images. 
None of the instruments can give a satisfactory image of "N" 
: the resolution of NRH is too low, and GMRT, in spite of its 
high resolution, does not image the bright core of ' W" properly. 
The resolution of this composite image is 49" as for the other 
images. The nns dynamic range of this image is 280, and the 
max dynamic range is 14. 

Upon an examination of high time resolution NRH data, it 
is evident that the relatively faint halo near 'W on the compos- 
ite images is in fact a group of short, bright and almost unpo- 
larized type HI bursts, beside the noise storms "S " and "N". In 
addition to the short time scale motions of "S " and "A^' men- 
tioned earlier, this shows that the observations presented here 
are only illustrative of the possibilities of the joint instrument. 
Higher time resolution observations with this joint instrument 
are very much needed in order to derive physical interpretations 
of solar radio emissions in this frequency range. 



5. Summary and conclusion 

We have demonstrated the capablities of a new NRH + GMRT 
composite instrument for obtaining snapshot images of the sun 
at 327 MHz. We have shown examples from simulations with 
two models (of differing complexity) of the radio sun, and from 
real observations of complex, evolving noise storm sources on 
Aug 27, 2002. 

The simulations show the capabibility of the joint instru- 
ment, to produce snapshot images of a complex sun, with struc- 
ture sizes ranging the resolution (a 50 ) up to the size of the 
whole sun (ss 40 ). The main conclusions from the simulations 
are : 

- GMRT and NRH are highly complementary : long GMRT 
baselines provide resolution, whereas dense NRH uv cov- 
erage at short baselines prevent aliasing of large structures 
which can occur with GMRT alone. 

- The composite instrument allows a substantial advance in 
snapshot imaging of complex objects. The resolution and 
dynamic ranges of composite images are far better than 
those of images from the individual instruments (Fig 7). 
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- The quality of the calibration can limit the resolution and 
produce appreciable artefacts. 

We produced composite 17 sec snapshot images (from ac- 
tual observations of the sun on Aug 27 2002) of structures be- 
tween 60 " and ~ 200 " in size with a resolution of 49 " and 
rms dynamic ranges of 250-420. The quality of the compos- 
ite image is far better than those of images from the individual 
instruments (eg fig. 8). The max dynamic ranges of the snap- 
shot maps are > 20. To the best of our knowledge, these are 
the highest dynamic range snapshot maps of the sun at meter 
wavelengths. Until now, high dynamic range radio maps of the 
sun were typically made by synthesis imaging over time peri- 
ods of a few hours. As mentioned in the introduction, high dy- 
namic range images would be essential in studying phenomena 
like bright radio bursts occurring along with (fainter phenom- 
ena hke) coronal mass ejections. 

These observational results are only illustrative of the ca- 
pabihties of the technique we have employed. Merging GMRT 
and NRH data will potentially yield ^ 1000 basehnes per snap- 
shot and would yield very high resolution, high dynamic range 
maps. This is roughly equivalent to simultaneously using the 
capabihties of the VLA B, C and D arrays (the comparison con- 
cerns only the number of basehnes, since these VLA configu- 
rations cannot be used simultaneously to produce snapshots). 

The resolution of the maps presented in this paper was lim- 
ited to 49 ", which is well below what one would crudely expect 
from the GMRT. This hmitation in the resolution is mostly be- 
cause: 

- We are dealing with snapshot images and the uv coverage 
of the GMRT is sparse, especially for large baselines. The 
deconvolution procedure does not work well for complex 
objects under these conditions. We need to introduce sub- 
stantial tapering of the visibilities (see § 3.2.2) in order to 
limit artefacts. This is a basic limitation on snapshot im- 
ages, which would not concern synthesis composite im- 
ages because the uv coverage (for synthesis images) would 
be much less sparse, especially for large baselines. When 
the object is simple enough (eg fig. 9), however, one can 
achieve (from snapshot data) the full resolution correspond- 
ing to the size of GMRT by using classical model fitting. 

- We had to flag several GMRT baselines owing to malfunc- 
tioning antennas and radio frequency interference. 

- We encountered calibration problems with GMRT data. We 
have checked the eff'ect of cahbration errors with simula- 
tions, and it is clear that they make a significant difference 
in the level of clean artefacts in the final images. The res- 
olution of 49 we have used is the best compromise be- 
tween achieving the best possible resolution and maintain- 
ing a reasonably low level of clean artefacts. 

For future observations, it is desirable to: 

- Achieve the best possible calibration, especially for the 
GMRT. The GMRT is probably more sensitive than NRH 
to phase deviations because of its sparser uv coverage. 

- use the shortest possible integration time for GMRT obser- 
vations (ss 2 sec, rather than 17 sec. as in the current obser- 
vations). This is because radio bursts can vary on timescales 
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much shorter than 17 sec, and the data presented here could 
have integrated over several potentially interesting time- 
varying phenomena (as was the case for the snapshot at 
09:04:04, Aug 27 2002). 

This could lead to images with better resolution and dy- 
namic ranges. We could then conclusively answer some long- 
standing questions in solar physics, as mentioned in the intro- 
duction. 
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